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Aggarwal  et  al.  [A.  Aggarwal,  M.M.  Klawe,  S.  Moran,  P.  Shor,  R.  Wilber, 
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every  vertex  of  a  convex  n-gon.‘  This  thesis  extends  the  results  of  Aggarwal  et.  al. 
by  developing  the  following  algorithms: 

l)  An  optimal  algorithm  to  find  all  farthest  vertices  for  every  vertex  of  a  convex 
polygon^ 

2^  An  0(kn  log  k)  time  algorithm  to  find  k  farthest  vertices  for  every  vertex  of  a 
convex  n-gonj  f  5?  r, 

3)  An  0(r&)  algorithm  to  sort  the  distances  of  all  the  vertices  of  a  convex  n-gon 
with  respect  to  each  vertex  of  the  convex  n-gon] 

4}  A  worst-case  optimal  algorithm  to  sort  a  set  of  numbers  given  lower  bounds 
on  the  ranks.'N 
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Introduction 


Given  a  set  5  of  n  points  in  the  plane,  for  each  p,  €  5  we  want  to  find  a  point 
p.  6  5.  j  it  >,  such  that 

lr 

d(.P,.  Pj)  =  max  d(p„  pi). 

0</<n-l 

where  d(pi,  pj)  is  the  Euclidean  distance  between  p,  and  p}.  We  will  cadi  this 
the  farthest  neighbor  problem.  Toussaint  and  Bhattacharya  [TB81],  and  Preparata 
'P771  have  shown  a  0(n  log  n)  time  bound  for  this  problem.  However,  an  Q(  n  log  n  j 
time  bound  is  obtained  under  the  assumption  that  the  input  to  the  algorithm  is  an 
arbitrary  set  of  points  in  the  plane.  If  we  restrict  the  input  to  be  the  vertices  of  a 
convex  polygon,  then  Aggarwal  et.  al.  [AKMSWS7]  have  shown  a0(n)  time  bound. 
In  section  2  we  define  monotone  matrices  that  are  used  in  [AKMSW87]  to  solve  the 
farthest  neighbor  problem.  We  use  these  matrices  to  obtain  our  results  as  well. 

We  extend  their  results  in  the  following  ways: 

1.  The  all  farthest  neighbors  problem  for  a  set  S  of  n  points  is  to  find  for  every 
point  pi  €  S,  all  the  points  Pj,  j  ^  t,  such  that 

d(Pi,  Pj)  =  0<max  pi). 

In  section  3.1  we  show  how  to  find  all  farthest  neighbors  for  every  vertex  of 
a  convex  polygon  and  in  section  3.2  we  prove  0(n)  time  complexity  for  our 
algorithm.  In  section  3.3  we  give  an  application  of  our  algorithm. 

2.  The  k  farthest  neighbors  problem  for  a  set  5  of  n  points  is  to  find,  for  every 
point  p,  6  S,  the  set  of  points  Pi,  |P,j  =  k,  such  that  every  pj  6  P,  is  farther 
from  pi  than  any  pi  £  P,.  In  section  4.1  we  show  how  to  find  k  farthest 
neighbors  for  every  vertex  of  a  convex  polygon  and  in  section  4.2  we  prove 
O(knlogk)  time  complexity  for  our  algorithm.  Section  4.3  contains  some 
comments  on  the  lower  bound  for  this  problem. 

3.  The  sorting  problem  for  a  planar  point  set  is  as  follows:  given  two  sets  S  and 
T  containing  n  and  m  planar  points,  respectively,  sort  the  points  of  T  with 
respect  to  their  distance  from  p,  for  each  p  6  5.  If  T  =  5.  then  the  problem 
reduces  to  sorting  5  with  respect  to  itself. 

Clearly,  we  can  sort  T  with  respect  to  S  in  O(nmlgm)  time,  using  standard 
sorting  techniques1.  This  naive  approachT however,  is  not  optimal  when  5  is 
a  convex  set,  i.e.  it  contains  the  vertices  of  a  convex  polygon.  In  Section  5.1. 
we  show  how  to  sort  T  with  respect  to  5  if  5  is  a  convex  set.  In  section  5.2 
we  prove  0(nm  4-  m2)  time  complexity  for  our  algorithm,  which  implies  that 

'[CIR89],  put  II. 
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we  can  sort  5  with  respect  to  itself  in  0(n2)  time.  Sorting  a  convex  set  S 
with  respect  to  itself  corresponds  to  the  case  k  =  n  in  the  k  farthest  neighbors 
problem  if  we  require  that  the  list  of  farthest  neighbors  be  produced  in  sorted 
order. 

Section  5.3  contains  some  comments  on  the  lower  bound  for  this  problem. 

We  -extend  this  result  in  Section  5.4.  obtaining  an  O(n2logf)  time  algorithm 
for  sorting  an  arbitrary  point  set  S  (with  l  convex  layers)  with  respect  to  itself. 

In  addition  to  the  problems  concerning  finding  farthest  neighbors,  in  section 
6  we  show  how  to  sort  a  set  of  numbers  gi  ven  a  lower  bound  on  the  tank  of -each 
number  in  the-sorted  order.  We  give  an  algorithm  that  solves  this  prohlem  in  section 
6.2.  find  the  time  complexity  of  the  algorithm  in  section  6.3,  .and  prove  that  it  is 
optimal  in  the  worst  case  in  section  6.4.  This  problem  is  related  to  other  questions 
of  sorting  using  partial  information  which  are  investigated  in,  for  example,  [BT80], 
[F76],  [GMPR77],  [KS84],  (LS83J. 

All  the  lower  bounds  in  this  paper  are  given  in  the  algebraic  decision  tree  model  of 
computation2.  Briefly,  in  this  model  the  cost  of  running  the  algorithm  is  considered 
to  be  proportional  to  the  number  of  comparisons  along  the  longest  path  from  the 
root  .to  a  leaf  in  the  decision  tree  that  corresponds  to  the  execution  of  the  algorithm. 


’For  detailed  discussion,  tee  page  30  of  [PS85]. 
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2  Monotone  Matrices 

In  sec*, on  we  explain  the  construction  used  in  [AKMSW87;  and  by  us  to  solve 

'lie  farthest  neighbor  problems.  A  two-dimensional  matrix  M  =  is  called 

monotone  if  "tie  maximum  value  in  the  i-th  row  lies  below  or  to  the  right  of  the 
maximum  value  in  the  u  -  l)-st  row.  If  a  row  has” several  maxima  then  we  will  take 
the  leftmost  one.  A  matrix  M  is  called  totally  monotone  if  every  2  x  2  submatrix 
( i.e. .  every  2x2  minor)  is  monotone  (See  Figure  l). 


Figure  1 :  Every  2x2  minor .  given  by  a,b,c,  and  d,  of  a  totally  monotone  matrix 
is  monotone,  i.e.  it  is  not  possible  that  a<b  and  c>d. 

Aitnougn  the  question  of  finding  the  row  maxima  in  a  two-dimensional  totally 
monotone  matrix  may  seem  rather  odd  at  first  glance.  [AKMSW87]  and  [AP881 
show  that  a  wide  variety  of  problems  can  be  reduced  to  one  or  more  instances 
of  this  problem.  For  example,  the  totally  monotone  property  arises  in  problems 
that  deal  with  polygons  whose  vertices  obey  the  quadrangle  inequality.  Suppose 
a  convex  polygon  has  n  vertices  numbered  0, 1..., n  —  1.  If  we  take  any  4  distinct 
vertices,  »'j,  « j,  ji,  and  jj,  such  that  0  <  »'t  <  ij  <  j\  <  ;i  <  «  -  1,  and  form  a 
quadrangle  then  by  the  quadrangle  inequality,  the  sum  of  lengths  of  the 

diagonals  d(i\,ji)  +  d(ij,ji)  is  strictly  greater  than  the  sum  of  the  lengths  of  the 
sides  d(»i,;3)  •+-  d(i?,ji)  (See  Figure  2). 


h  h 


a 

b 

c 

T 

We  define  a  matrix  .V/  =  (m,., }  corresponding  to  a  convex  polygon  as  follows 


.  See  Figure  3):  Vi.  0  <  i  <  L.  m,  :  =  d(i.<j  mod  n )  j  for  i  <  j  <  <  -  n.  and 
m  =  -  x  for  all  ot  her  j . 
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Claim  1  Matrix  M  13  a  totally  monotone  matrix. 

Proof  Take  any  2x2  minor  (See  Figure  l).  If  none  of  the  entries  in  the  minor  are 
-x  and  »i  <  >2  <  J\  <  ;'j.  then  it  is  not  possible  that  a  <  b  and  c  >  d  because  of 
the  quadrangle  inequality.  »i  <  i 3  and  <  jj  have  to  hold  simply  from  the  way  we 
have  drawn  the  minor.  If  >  ji,  then  c  =  -x.  Similarly,  if  any  other  inequality 
among  »i  <  17  <  J\  <  h  does  not  hold,  then  at  least  one  of  the  entries  in  the  minor 
is  -x.  which  is  the  case  we  consider  next. 

If  b  js  -x  or  c  as  -x,  then  either  a  <  6  or  c  ?  d.  If  a  =  - X .  then  by  the 
way  we  have  defined  A/,  either  b  =  -xore  =  -x  and  the  previous  case  applies. 
Similarly,  if  d  s  -x,  then  either  c  =  - x  or  b  *  -x.  Thus,  M  is  monotone.  ■ 

If  each  entry  of  a  totally  monotone  nxm  matrix  M  can  be  computed  in  constant 
time,  then  Aggarwal  et  al.  [AKMSW87]  showed  how  to  find  the  leftmost  maximum 
in  each  row  of  an  n  x  m  matrix  A/  in  0(n  4-  m)  time  on  a  sequential  RAM.  In 
the  following  discussion  we  will  refer  to  the  [AKMSW87]  algorithm  as  the  SMAWK 
algorithm,  following  the  convention  used  by  the  authors.  .Vote  that  the  SMAWK 
algorithm  does  not  explicitly  create  the  entire  matrix  XI  (that  would  take  0(nm ) 
time):  rather,  it  only  computes  0(n  +  m)  entries  of  M.  Thus,  the  [AKMSW87] 
algorithm  can  find  the  leftmost  farthest  neighbor  of  every  vertex  on  a  convex  n-gon 
in  O(n)  time. 

Before  concluding  this  section,  we  note  that  if  /,  denotes  the  number  of  farthest 
neighbors  of  the  i-th  vertex,  then  it  can  be  shown  (see  appendix)  that  £*J0l  /,  <  2n. 


Observe  that  there  are  convex  polygons  for  which  />  -  2n.  For  example,  put 
n  -  1  vertices  on  an  arc  such  that  there  are  vertices  at  both  endpoints  of  the  arc  and 
the  sire  of  the  arc  is  |  of  a  circle.  Let  the  center  of  the  circle  be  the  n-th  vertex. 
The  center  has  n  -  1  farthest  neighbors,  each  endpoint  of  the  arc  has  2  farthest 
neighbors,  and  each  of  the  other  n  -  3  points  on  the  arc  has  1  farthest  neighbor. 
The  sum  is(n  -  l)  +  2x2  +  (n-3)x  1  =  2 n. 

Two  vertices  »  and  j  are  called  symmetric  farthest  neighbors  if  f(i)  =  j  and 
f(j)  —  »•  Since  f,  <  2 n.  there  are  at  most  n  symmetric  farthest  neighbor 

pairs.  The  polygon  described  above  again  achieves  the  upper  bound  of  n. 


3  The  All  Farthest  Neighbors  Problem  for  a  Convex 
Polygon 


3.1  The  Algorithm  for  Finding  All  Farthest  Neighbors  of  a  Convex 
Polygon 

Let  .V/  =  be  4  totally  monotone  n  x  2n  matrix  corresponding  to  a  convex 

polygon  that  was  defined  in  section  2.  We  use  the  SMAWK  algorithm  to  find  the 
leftmost  maximum  in  each  row  of  .V/,  Now.  consider  a  different  totally  monotone 
matrix  .V/'  such  that  Vi.  0  <  i  <  n- 1.  m\ }  =  d((n  -  1  -i)  mod  n.(2n- 1  )  mod  n) 

for  i  <  ;  <  i  f  n  and  m't  J  =  -  x  for  all  other  j. 

Note  that  the  matrix  A/'  is  simply  \i  rotated  upside  down,  or  equivalently,  M 
flipped  horizontally  and  vertically  (See  Figure  4). 


Figure  4.  The  matrix  from  Figure  1  flipped  horizontally  and  vertically. 


Such  a  rotation  preserves  the  totally  monotone  property:  it  is  not  possible  to 
have  d  <  c  and  b  >  a  (which  is  the  same  as  saying  that  it  is  not  possible  to  have  a  <  b 
and  c  >  d).  Again,  we  use  the  SMAWK  algorithm  to  find  the  leftmost  maximum  in 
each  row  of  A/'.  But,  by  the  way  we  have  defined  A/\  the  leftmost  maximum  in  row 
«  of  .V/'  is  the  rightmost  maximum  in  row  (n  -  1  -  i)  mod  n  of  A/.  Therefore,  we 
now  have  the  rightmost  maximum  in  each  row  of  A/.  Last,  for  each  row  i,  check  all 
the  entries  between  the  leftmost  and  the  rightmost  maxima  in  i  for  other  maxima 
in  that  row.  The  summary  of  the  algorithm  follows. 


This  section  is  joint  work  with  Aiok  Aggatwal. 
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run  SMAWK  on  M 
for  i  —  l  to  n  do 

/,  —  column  of  the  leftmost  maximum  ...  .^.v  i  of  M 
output  (i./, ) 
run  SMAWK  on  M' 
for  i  —  1  to  n  do 

r,  —  column  of  the  leftmost  maximum  in  row  ( n  -  1  -  t)  mod  n  of  M' 
if  /,  jt  r;  than  output  i.i.ri) 
for  t  —  1  to  n  do 

for  ;  -  +  1  to  r,  -  1  do 

if  m,  j  =  m, then  output  (i.j) 

3.2  Time  Complexity  of  the  Algorithm 

The  first  two  loops  of  the  algorithm  and  the  two  executions  of  the  SMAWK  algorithm 
take  O(n)  time.  Clearly,  all  the  row  maxima  have  to  be  between  the  leftmost  and 
the  rightmost  maxima  in  that  row.  Note  that  m,  /t  =  m,  r, .  All  we  need  to  prove 
now  is  that  the  third  loop  of  the  algorithm  takes  0(n)  time. 

From  the  quadrangle  inequality  we  know  that  a  +  d  >  6  +  c  (See  Figure  2).  From 
this  we  derive  the  following  observations: 

a  <  b  implies  c  <  d 
c  >  d  implies  a  >  b. 

Lemma  1  r,  <  min{2n  -  l./.+j}. 

Proof  Suppose  r,  >  mm{2n  -  l,i,+i}.  Let  z  -  m,> ,  =  to,-./,,  y  =  m, ./,*,, 

*  =  and  w  =  (See  Figure  5). 


Figure  5  ~ 

Since  x  is  the  maximum  value  for  row  s,  it  must  be  true  that  y  <  x.  Then,  by 
our  first  observation,  z  <  w,  which  contradicts  the  fact  that  z  is  the  maximum  for 
row  i  +  1.  Thus,  r,  <  m«n{2n  -  1  ■ 
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Lemma  2  The  third  loop  of  the -algorithm  takes  0(n)  time. 

Proof  Let  Tv,  be  the  region  between  mlJt  and  m,.r,.  By  Lemma  1.  Vi,j.  i  j. 
Tv,  and  TZ,  cannot  overlap  horizontally.  Thus,  each  x-coordinate  of  matrix  M  will 
be  checked  in  at  most  one  such  region  in  the  third  loop  of  the  algorithm.  Since  the 
horizontal  dimension  of  M  is  2n,  we  need  to  check  at  most  0(n)  entries  of  M.  W 

Since  we  output  at  least  one  farthest  neighbor  for  each  of  the  n  vertices,  the 
O(n)  running  time  of  our  algorithm  is  optimal. 

.3.3  An  Application:  All  Symmetric  Farthest  Neighbors  of  a  Simple 
Polygon 

It  is  well  known  that  the  farthest  neighbor  of  any  point  inside  a  convex  region  is 
one  of  the  vertices  on  the  convex  hull.  Also,  the  convex  hull  of  a  simple  polygon 
can  be  found  in  O(n)  time  ([GY83]).  So,  to  find  all  symmetric  farthest  neighbors  of 
a  simple  polygon  we  first  find  the  convex  hull  of  the  simple  polygon.  Then,  using 
the  algorithm  in  section  3.1,  we  find  all  farthest  neighbors  of  the  convex  hull,  and 
finally  scan  the  list  of  vertices  ofihe  convex  hull  to  find  which  ones  are  symmetric 
farthest  neighbors.  Thus,  all  symmetric  farthest  neighbors  of  a  simple  polygon  can 
be  found  in  @(n)  time.  This  settles  a  problem  raised  in  [T83]. 
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4  The  k  Farthest  Neighbors  Problem  for  a  Convex 
Polygon 

4.1  The  Algorithm  for  Finding  k  Farthest  Neighbors  of  a  Convex 
Polygon 

Let  M  =  {m,j}  be  a  totally  monotone  n  x  m  matrix.  Our  algorithm  consists, 
roughly  speaking,  of  three  parts.  In  the  first  and  second  parts,  we  decompose  each 
row  of  Af  as  follows:  for  each  row  j,  we  partially  order  elements  of  row  j  into  2 k 
sorted  lists  each  with  upto  k  elements  such  that  the  k  largest  elements  of  row  j  are 
guaranteed  to  be  among  the  lists  for  row  j1.  In  the  third  part,  for  each  row  j  we 
find  the  k  largest  elements  in  j.  We  achieve  this  by  partially  merging  the  2k  lists 
for  j  that  we  found  in  the  first  two  parts. 

Let  e[i,j]  be  the  column  of  the  largest  element  in  the  s'-th  sorted  list  for  row  j. 
The  first  part  has  k  phases.  In  the  s-th  phase  we  find  (2i  -  l)-st  and  2i-th  list  of 
upto  k  elements  for  each  row  j  using  the  2(i  -  1)  lists  found  in  the  previous  phases. 
We  achieve  this  by  running  the  SMAWK  algorithm  on  the  matrix  M',  where  M'  is 
M  with  columns  c[i',/]  for  1  <  «'  <  i  -  1  and  1  <  j'  <  n  removed.  Running  the 
SMAWK  algorithm  gives  us  c{$,j]  for  each  row  j.  Moreover,  because  M  and  M' 
are  totally  monotone,  c[i,j]’s  actually  give  us  more  information.  If  we  look  at  the 
matrix  At  given  by  the  columns  c[»,  1]  through  c{i,  n],3  we  observe  that  the  largest 
element  in  row  j  of  M  is  Furthermore,  each  row  of  M  is  bi tonic. 

Lemma  3  Vj  and  Vi, 

mi.eCi,U  S  •' ■  <  a™* 

mA4Ml  -  m;.e(M+ll  -  ’ "  - 

Proof  Suppose  not,  i.e.,  there  is  some  i  and  some  row  j  of  Af  such  that  for  some 
a  and  t,  t  <  s  <  j,  my^j  >  my^j  (the  case  of  j  <  »  <  t  follows  analogously).  If 
c[i,a]  *  c(M],  then  certainly  =  my^.  If  c[»,s]  jt  c(«,t],  then  since  m,ieM 

is  the  largest  dement  in  row  s,  <  m^j.  But  if  Af  is  totally  monotone,  then 

it  is  impossible  that  and  my^j  >  my^j.  Thus,  our  assumption 

most  have  been  wrong.  ■ 

Thus,  for  each  row  j,  the  c(i,y]’s  define  2 k  sorted  lists: 

TUs  ssctioa  is  joist  work  with  James  Park. 

‘Not*  that  ws  seem  to  be  deaKag  with  Ofk* »)  alemaata  is  these  lasts.  We  oaly  see  these  lists 
for  the  explaaatioe;  as  we  shall  see,  the  algorithm  actaally  ases  oaly  O(kn)  elemeato  of  if. 

3 Note  that  c{»,  j]  may  aqeal  c{»,  j  +  1)  for  some  >’»,  bet  each  colama  of  if  appear*  oaly  oace  ia 

St. 
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^  J-i)  2:  • •  •  > 
mj,c{l  jJ  -  m>,c{lJ+l]  -  "  '  - 

1  -  ^  '  ^  mi.<42.i] 

-  mi.c(2j+l]  2  "  >  mj43.m] 

mj,e(M  -  ^  ‘  •  2  mj.efJU] 

mi4*j]  ^  mi4*o+M  -  - 

Unfortunately,  we  ran  into  a  problem  here  since  c[i,  j']  may  be  equal  to  c[t,/  + 1] 
for  some  number  of  (j')a.  This  is  a  problem  because  in  order  to  achieve  the  desired 
time  bounds  we  want  each  member  of  each  list  to  be  a  different  entry  of  M .  To 
solve  the  difficulty  we  need  to  “skip”  the  repeating  entries.  This  is  accomplished  in 
the  second  part  by  defining  the  skip  variables 

=  j-  I  if  :  1  <  ]'  <  j  and  «[*»/]  =  c[i,/]}|  -  1 

dovm[iJ]  -j  +  |{/  :  n>j’>  j  and  c[i,/]  =  c[t,y]}|  +  1 
Now,  we  can  define  2k  sorted  lists  for  each  row  j  in  which  no  entry  of  M  is  repeated: 

-  *  *  * 

mi.c(2j]  -  mi.c(2.«sl2jD  -  ^  ■ 

^  m>,c( 2,doum[2jjj>  ”»i42.<<o«m[2,«iown(2j^  *  *  * 


Also,  observe  that  Vj,  m^ij]  >  mj^j]  >  In  other  words, 

“♦’dd  be  the  i-th  largest  in  row  j  only  if  is  one  of  the  t  -  1 

largest. 

In  the  third  part  of  our  algorithm,  for  each  row  j  of  Af  we  merge  parts  of  these 
lists  to  get  the  k  largest  elements.  This  part  of  the  algorithm  consists  of  k  phases 
(for  each  row  j).  In  phase  i  we  find  i-th  largest  value  in  row  j  using  i  -  1  largest 
values  in  that  row.  We  partially  merge  these  lists  using  a  data  structure  H.  H  must 
support  the  following  operations:  insert,  find  the  maximum  element  and  delete  the 
maximum  element.  One  suitable  implementation  for  is  a  heap9.  The  algorithm 
is  summarized  below. 


’For  detailed  diacunon,  Me  chapter  7  of  [CLRS9]. 
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Af'«-  M 

for  i  *-  1  to  k  do 

run  SMAWK  on  M • 
for  j  *-  1  to  n  do 

c[»,j]  <—  column  of  the  maximum  in  row  j  of  Af' 
for  j  1  to  n  do  ” 

M'  *-  M'  -  column  c[»,  j] 
for  i  <-  1  to  k  do 

down  «-  1;  up «-  1;  high  *-  1 
while  high  <  n  do 
low «—  high 

while  e[»,  high]  =  c[»,  high  +  1]  and  high  <  n  do 
high  <—  high  4- 1 
for  j  low  to  high  do 
up{i,  ji]  «-  low  —  1 
doum[i,  j]  ♦-  high  +  1 
high  4-  high  +  1 
for  j  *-  1  to  n  do 
L[\,j]  4- 
initialize  H 

insert(F,  mi(e(li<town{W]],  roi>c( 2jj) 

for  i  4-  2  to  Jb  do 

mi4»'j'l  *"  ftnd-max(iJ) 

delete- max(  if) 

if  j'  <  j  then  insert(£T, 

if  j'  >  j  then  insert  (if, 

if  /  =  3  then  insert  (if,  miiC(i,+1>i]) 

4.2  Time  Complexity  of  the  Algorithm 

The  first  part  of  the  algorithm  is  running  the  SMAWK  algorithm  k  times  and  thus 
takes  0(fc(n  +  m))  time.  The  second  part  does  a  constant  number  of  operations 
for  each  e(i,j],  thus  taking  0(k(n  4-  m))  time.  In  the  third  loop,  for  each  i,  we 
are  running  k  insert,  find-max  and  delete-max  operations  on  a  heap  H.  Each  of 
these  operations  takes  no  more  than  0(log|Jf|)  time,  where  |27|  is  the  number  of 
dements  in  if  at  the  time  the  operation  is  performed.  Luckily,  H  never  has  too 
many  dements.  In  fact,  |lf|  <  2k  since  at  any  point  H  contains  at  most  one  dement 
from  each  of  the  2k  lists.  Thus,  the  second  part  of  the  algorithm  takes  0(kn  log  k) 
time  and  the  total  time  for  the  algorithm  is  0(kn  log  k). 

Note  that  this  algorithm  could  be  used  to  find  k  farthest  neighbors  on  a  convex 
n-gon  in  O(knlogk)  time. 


13 


4.3  Some  Comments  on  the  Lower  Bounds 

To  output  n  lists  of  k  sorted  elements  each  would  seem  to  require  fl(kn  logib)  time. 
However,  these  lists  are  not  necessarily  independent.  In  terms  of  our  monotone 
matrix  M ,  these  lists  are  not  necessarily  column  independent,  i.e.,  one  of  the  k 
largest  elements  in  row  t  may  be  in  the  same  column  as  one  of  the  k  largest  elements 
in  row  i'.  In  terms  of  a  polygon,  the  same  vertex  could  be  included  as  one  of  the 
k  farthest  vertices  for  more  than  one  vertex.  Yet,  no  o(kn logib)  time  algorithm  is 
known  that  gives  the  k  largest  elements  even  when  the  elements  are  not  required  in 
sorted  order. 
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The  Sorting  Problem  for  a  Planar  Point  Set 


5.1  The  Algorithm  for  Sorting  the  PQints  of  an  Arbitrary  Set  T 
with  respect  to  a  Convex  Set  S 

Let  a  convex  set  S  be  given  by  the  vertices  of  a  convex  polygon  Ps-  Let  {a®,  *  •  • » ^n-i } 
be  the  points  of  5  in  clockwise  order  around  the  perimeter  of  Ps  and  let  {to,  •  •  • ,  tm-i } 
be  the  points  of  T  in  any  order.  Denote  by  m,  j  the  distance  from  a<  to  t,  and  by 
j[it  r]  the  index  of  the  point  in  T  whose  distance  from  a,  is  r-th  smallest  among 
points  in  T.  (For  convenience,  we  define  j[0,  r]  =  r  for  1  <  r  <  m.) 

Our  algorithm  consists  of  n  phases,  where  in  the  t-th  phase,  we  compute  j[i,r] 
for  1  <  r  <  m,  using  the  values  j{»  —  1,  r]  we  computed  in  the  previous  phase.  We 
do  this  using  an  insertion  sort,  inserting  then  *>[i_i,2]>  then  *;(,•_  1,3] >  and  so 

on  through  This  is  summarized  below. 


for  r  «-  1  to  m  do 
j[0,r]  «-  r 

for  *  0  to  n  -  1  do 

for  r  «-  1  to  m  do 

iM  1*rl 

k «-  r 

while  k  >  1  and  m,^  >  do 

*-  j[i,k  -  1} 
j[i,k-l]*-t 
k  «—  k  —  1 


5.2  Time  Complexity  of  the  Algorithm 


We  claim  that  the  i-th  phase  of  our  algorithm  requires  0(m+li)  time,  where  /,-  is  the 
number  of  inversions  that  occur  in  going  from  the  ordering  of  T’ s  points  by  distance 
from  s,_i  to  the  ordering  of  T»  points  by  distance  from  a,-.  An  inversion  corresponds 
to  a  pair  of  points  (tj,ty)  from  T  such  that  m,_jj  <  but  m,j  >  m,  j>.  Since 

the  perpendicular  bisector  of  any  pair  of  points  (t>,  ty)  from  T  intersects  Ps  at  most 
twice, 


II— 1 


£/<  <  2 


Thus,  our  algorithm  runs  in  0(nm  +  m3)  total-  time. 
This  aactioa  is  joist  work  with  Jams  Park. 
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Note  that  this  algorithm  can  be  used  to  sort  the  rows  of  an  n  x  m  totally 
monotone  matrix  M  =  {m,j}  in ~0(nm  +  m2)  time,  as  total  monotonicity  implies 
that  for  each  pair  of  columns  j  <  j ',  either 

(1)  Vi  m,  j  <  m,v,»  or 

(2)  3t,  0  <  t  <  n  -  1  such  that  m,<j  >  mi>j>  for  i'  <  i  and  <  m,»j*  for  *'  >  t. 

Thus,  we  can  sort  the  set  S  of  vertices  of  a  convex  polygon  with  respect  to  itself 
in  0(n2)  time. 

5.3  Some  Comments  on  Lower  Bounds 

Sorting  T  with  respect  to  5  would  seem  to  require  0(nm)  time.  This  is  because 
specifying  n  orderings  of  m  points  in  the  obvious  way,  namely,  by  outputting  n 
ordered  lists  of  length  m,  itself  requires  n(nm)  time.  However,  this  is  not  necessarily 
the  case.  Certainly,  we  can  reduce  the  problem  of  sorting  m  numbers  to  the  problem 
of  sorting  m  paints  of  T  with  respect  to  one  point  of  S.  To  do  this  reduction,  we 
take  the  set  of  m  numbers  {zo<**'.*m-i}  and  we  create  points  cd 

T  as  follows:  t,-  =  (x,-,0)  for  0  <  *  <  m  —  1.  Furthermore,  we  let  *o  =  (0,0). 
Then,  sorting  T  with  respect  to  so  is  equivalent  to  sorting  {zo.’*'.Xm-i}-  Thus, 
we  know  that  sorting  T  with  respect  to  S  requires  O(mlogm)  time.  However,  once 
we  have  sorted  T  with  respect  to  so.  it  is  not  clear  if  sorting  T  with  respect  to  si 
requires  as  much  time.  The  difficulty  with  showing  the  ft(nm  +  m2)  lower  bound 
for  sorting  T  with  respect  to  5  is  that  only  0(n*)  differing  orderings  of  T’s  points 
are  possible  —  the  (,)  perpendicular  bisectors  of  pairs  of  pants  from  5  divide  the 
plane  into  0(n4)  regions,  and  there  is  a  one-to-one  correspondence  between  regions 
and  orderings.  Thus,  it  may  be  possible  to  specify  an  ordering  of  T’s  points  with 
respect  to  some  sj’e  in  o(n)  time,  which  suggests  that  a  o(nm)  time  solution  to  the 
problem  of  sorting  T  with  respect  to  5  is  still  attainable. 

5.4  An  0(n2  lg  £)  Time  Algorithm  for  Arbitrary  S 

We  can  solve  the  sorting  problem  for  an  arbitrary  set  S  of  n  planar  points  using  our 
0(nm  +  m2)  time  algorithm  for  sorting  an  arbitrary  set  with  respect  to  a  convex 
set.  We  use  two  partitions  of  $:  we  partition  5  into  l  subsets  Si,..., St,  each 
corresponding  to  a  convex  layer  of  5,1  and  l  subsets  5} , . . . ,  S},  each  of  size  n/l.  For 
1  <  *  <  l  and  1  <  j  <  l,  we  sort  Sj  with  respect  to  5,-  in  0(n,n/f  +  n2//2)  time, 
where  n,-  is  the  size  of  5,.  The  total  time  required  is 

ssoH+!)=so(-n+T)=oM- 

'We  cu  Sad  convex  layer*  of  n  point*  ia  0(nlog»)  time  ([PS85],  p.  16S). 
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For  each  point  p  of  S,  we  now  have  t  sorted  lists,  corresponding  to  the  points  of 
5{, . . .  ,Sf,  respectively.  As  we  can  merge  these  lists  in  0(n3 lgf)  total  time,  we 
obtain  an  0(n3  lgf)  time  algorithm  for  sorting  the  points  of  S. 


6  Sorting  Using  Lower  Bounds  on  the  Ranks 

6.1  Background 

Let  S  =  {so,  be  a  set  of  n  numbers.  To  sort  5  requires  ft(nlog  n)  time  and 

there  are  algorithms  that  achieve  this  lower  bound  ([CLR89],  part  II).  What  if  we 
are  given  some  partial  information  about  5?  The  lower  bound  may  not  necessarily 
hold.  Various  kinds  of  partial  information  about  S  have  been  researched  in  [BT80], 
[F76],  [GMPR77],  [KS84]  and  [LS83]. 

We  consider  the  following  problem.  Given  a  set  S  of  n  numbers,  we  have  the 
following  partial  information  about  the  final  sorted  order:  for  each  x,  €  S,  we  have 
a  lower  bound  /<  on  x,’s  rank,  i.eM  a  lower  bound  on  the  number  of  x,  €  S’  such  that 
Xj  <  x<.  How  much  time  does  it  take  to  sort  5? 

6.2  The  Algorithm  for  Sorting  5  Using  Lower  Bounds  on  the  Ranks 

We  create  n  buckets  (bj,  Bucket  bj  will  represent  rank  j  in  the  sorted  order 

of  5  and  by  the  end  of  the  algorithm,  6j  [0]  will  contain  the  y'-th  largest  element  of 
S.  We  put  into  bj  all  those  z,-’s  whose  lower  bound  on  the  rank  is  j .  Let  Cj  be  the 
number  of  x,-’s  placed  into  bj.  We  then  go  through  the  buckets  one  at  a  time,  from 
6i  to  bn.  We  use  a  data  structure  H  to  process  the  buckets.  To  process  bucket  bj, 
we  take  the  Cj  elements  in  bj  and  insert  them  into  H.  We  then  find  the  maximum 
element  in  H,  put  that  element  into  6,(0]  and  delete  it  from  H.  H  has  to  support  the 
following  operations:  insert,  find  maximum  element,  and  delete  maximum  element. 
As  in  section  4.1,  one  suitable  data  structure  for  H  is  a  heap.  The  algorithm  is 
summarized  below. 

for  j  «-  1  to  n  do 
Cj  -0 

for  t  <-  0  to  n  -  1  do 

ci,  cti  +  1 
-  *< 

initialize  H 

for  j  «-  1  to  n  do 

for  k  <-  1  to  Cj  do 
insert(27,  6/fifc]) 

Znu  *-  find-max(H) 

WO]  «“  Imu 

delete-max(H) 
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6.3  Time  Complexity  of  the  Algorithm 

Let  Q  =  CLi  (r,  -  /,•  +  1),  where  r,-  is  the  exact  rank  of  x,-  in  the  sorted  order2. 
Since  (ri,...,r„)  is  just  a  permutation  of  the  set  (l,...,n),  2Z"=1  r,-  =  l) .  Thus, 

Q  =  ^  -  E?=i  /.• 

Lemma  4  The  algorithm  takes  0{ n  log  £)  time. 

Proof  The  first  two  loops  take  0(n)  time.  Note  that  £"=1  <j-  =  n.  In  step  j  of 
the  third  loop,  we  insert  Cj  elements  into  H,  then  find  and  delete  the  maximum  in 
H.  Each  of  these  operations  takes  no  more  than  O(log  |/7j)  time,  where  \H\  is  the 
number  of  elements  in  H  at  the  time  the  operation  is  performed.  So,  if  we  let  w, 
be  the  number  of  elements  in  H  at  the  end  of  the  j-th  iteration,  then  during  the 
(j  +  l)-st  iteration,  we  pay  0(log(tOj)  +  log(u>j  +  1)  +  ...  +  log(u>j  +  ci+1  -  1))  to 
insert  Cj+i  elements  into  H  and  0(log(tOj  +  Cj+i))  to  find  and  delete  the  maximum. 
Note  that  our  definition  of  the  to/s  implies  that  wo  =  wn  =  0  and  for  0  <  j  <  n, 
0  <  Wj  <  n.  We  can  upper-bound  the  time  spent  on  the  (j  +  l)-st  iteration  by 
(«i+i  +  l)log(u?j  +  ei+1). 

Now  we  express  the  to/s  in  terms  of  the  c/s.  At  the  end  j-th  iteration  of  the 
third  loop,  H  contains  all  the  elements  we  have  inserted  during  the  first  j  iterations, 
minus  the  number  of  maxima  we  have  deleted.  Thus, 

i 

«i  =  £c.- j. 

iml 

Alternatively, 

i- 1 

=  (X>-0'-l)J  +  c,-l 

=  Wj.  I  +  Cj-  1. 

Then  Cj  =  Wj  -  wj.i  +  1.  So,  the  time  T  for  the  third  step  can  be  bounded  as 
follows: 


<  0  (±(cj  +  1)  log(u>j_1  +  Cj)  j 

^  0  +  1)  +  1)  +  (»f  -  t»i-i  +  1))  j 

<  0  +  2)le8(®>  +  l)j 


*Ne4s  that  we  an  oily  iMaf  r,’s  to  amplify  oar  explanation,  we  do  not  need  them  at  all  dace 
Q  is  amply  a  faactioa  of  b's  aid  a. 
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s  0  (•> lo*  *°  (Stt)) + + 0  (*—’■ l0*  (^+r)) 

+0(«nlog(w«  +  l))  +  0  ^J~jlog(u’j  + 1)^  * 

The  terms  corresponding  to  any  Wj  =  0  dissappear.  So,  henceforth,  we  assume  that 
Vj,  /  0.  Since  V;  Cj  >  0,  u>,_i  -  1  <  tn>,  and  thus,  ^  Then, 


r  <  o(»i^(:^))+^  (»>'<«  (5^))  +  --  +  0(»*-'lo«(^cr)) 

+0(wn  logfu^n  +  1 ))  +  0  +  l)j 

<  0  (mlog  (l  +  i))  +0  («log  (l  +  £))  +  ■••  +  o  («k-,l0*  (l  +  ^)) 


+O(ionlog(w»  + 1))  + 
Since  log(l  + 1)  <  0(()  and  v>n  *  0 


0  ^5Z1°8(u';  +  l)j  • 


T  <  o(isix^)+o(^x±)  +  ...  +  o(wn_1x^)+0^|;iog(«;i  +  l^ 
O  ^£log(«>  +  l)j  • 


<  0(n)  +  < 


Since  each  x<  remains  in  H  from  the  step  during  which  bucket  6|,  is  processed  until 
the  step  during  which  bucket  &r,  is  processed,  we  know  that  Y,%i  *>j  S  Q-  We  then 
want  to  find  the  maximum  value  of  » 

£log(i^  +  1)  *  log  +  1) 
j*  i  >*i 

subject  to 

£  wi  £  Q- 

The  is  achieved  when  »i  *  tpj  =:  ...  »  w*  =  {£.  Certainly  any  other 

solution,  in  particular  a  solution  that  obeys  the  constraint  ufo  *  0,  will  yield  a 
smaller  value.  Thus, 

Time  <  0(n)+  0  (lpg($  +  1)")  <0(nlog— ). 
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6.4  Worst-Case  Optimality  of  the  Algorithm 

Lemma  5  For  every  given  Q,  there  exists  a  sorting  problem  that  will  take  0(n  log  £) 
time. 


Proof  We  prove  by  contradiction.  Given  some.0,  suppose  that  there  exists  some 
algorithm,  A,  that  sorts  the  numbers  in  o(nlog»)  time.  Consider  the  following 
problem,  which  we  call  the  set-sorting  problem:  we  have  ^  independent  sets  of  ^ 
elements  each  and  we  want  to  sort  all  the  sets.  We  show  that  if  A  exists,  then  we 
can  solve  the  set-sorting  problem  in 


■(H-J) 


2 

time,  which  wonld  contradict  the  lower  bound  on  sorting:  to  sort  independent 
sets  of  £  elements  each,  requires 


n 


x  —  log 
n 


time. 

Our  algorithm  for  the  set-sorting  problem  works  as  follows.  First,  for  each  set 
Sj,  we  find  maxj,  the  maximum  element  in  Sj.  Then,  we  change  the  sets  so  that 
each  element  of  set  Sj  is  strictly  greater  than  each  element  in  set  Sj>  for  all  j'  <  j. 
We  accomplish  this  change  by  incrementing  each  element  of  Sj  by  maxj>. 
Observe  that  this  change  preserves  the  relative  order  of  elements  within  each  set  Sj. 
Now  we  assign  the  lower  bounds  on  ranks:  each  element  in  Sj  is  assigned  the  lower 
bound  of  +  1.  Thus,  we  have  formulated  the  set-sorting  problem  as  an  input 
to  A:  we  have  ^  x  £  =  n  elements  to  be  sorted  given  the  lower  bound  on  the  rank 
of  each  element.  Note  that  the  output  of  A  gives  us  a  solution  to  the  set-sorting 
problem.  More  specifically,  to  get  the  rank  of  any  element  in  Sj  we  take  the  rank  of 
this  element  as  given  by  A  and  we  subtract  The  summary  of  the  algorithm 

follows.  We  denote  the  i-th  element  in  Sj  by  x\j,  t]. 
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for  j  *-  1  to  ^  do 

rnaxj  < - oo 

for  i  «—  1  to  £  do 

if  x[j,i)  >  maxj  then  maxj  *-  *[;',*] 

incr  «—  0 

for  ; V  2  to  ^  do 

incr  <—  incr  +  maxj-i 
for  i  <—  1  to  £  do 

x[j,  i]  ♦-  x[j,  i]  +  incr 

run  A  on  *[y,*]’s  with  /(j,  i]*s 
for  j  «-  1  to  ^  do 

for  i  ♦—  1  to  £  do 


Now  we  want  to  analyze  the  running  time  of  this  algorithm.  All  the  loops  take 

0(^xf)-0W 

time.  Recall  the  definition  of  Q.  So,  for  algorithm  A, 

Qa  =  £  (rA\j,  •)  ~  l\j>  *]  +  !)• 


w 


We  claim  that 


Qa 


£  x  (£  +  l)  n7 

— -xr0^ 


Thus,  the  running  time  of  the  set-sorting  algorithm  is 

0(n)  +  o  ^nlog  ^  j  =  0(n )  +  o  (nlog  , 

which  is  the  contradiction  we  wanted  to  achieve. 


Thus,  onr  algorithm  is  optimal  in  the  sense  that  for  any  Q,  there  exists  a  sorting 
problem  that  reqnires  fi(nlog  £)  time. 
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7  Conclusion  and  Open  Problems 

We  have  shown  the  following  results  in  this  thesis: 

1.  An  optimal  algorithm  to  find  all  farthest  vertices  for  every  vertex  of  a  convex 
polygon. 

2.  An  O(knlogfc)  time  algorithm  to  find  k  farthest  vertices  for  every  vertex  of  a 
convex  n-gon. 

3.  An  0(n2)  algorithm  to  sort  the  distances  of  all  the  vertices  of  a  convex  n-gon 
with  respect  to  each  vertex  of  the  convex  n-gon. 

4.  A  worst-case  optimal  algorithm  to  sort  a  set  of  numbers  given  the  lower  bound 
on  the  rank  of  each  number  in  the  sorted  order. 

The  following  open  questions  remain: 

1.  A  Q(kn)  algorithm  to  find  k  farthest  vertices  for  every  vertex  of  a  convex 
n-gon. 

2.  A  proof  that  fl(n2)  is  the  lower  bound  for  sorting  the  distances  of  all  the 
vertices  of  a  convex  n-gon  with  respect  to  each  vertex  of  the  convex  n-gon. 

3.  An  information-theory  optimal  algorithm  to  sort  a  set  of  numbers  given  the 
lower  bound  on  the  rank  of  each  number  in  the  sorted  order. 
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Appendix 


Let  potPi,  —,Pn- 1  be  the  vertices  of  the  simple  convex  polygon  in  clockwise  order 
and  let  /,  be  the  number  of  farthest  neighbors  of  the  i-th  vertex. 

Theorem  1  /<  <  2». 

Proof  Observe  that  corresponding  to  each  point  p,  is  a  6,-,  0  <  b{  <  n  —  1,  and 
an  e,,  0  <  e<  <  n  -  1,  that  together  define  a  range  ft0P&i+i,— >ft,+e,  of  consecutive 
points  on  the  polygon  not  containing  p,-  (where  the  indices  are  taken  modulo  n)  such 
that 

(1)  pi,,  is  a  farthest  neighbors  of  p,-, 

(2)  ftj+e<  is  a  farthest  neighbor  of  pt,  and 

(3)  all  other  farthest  neighbors  of  p,  (if  any)  are  among  Pb,+i,-  •  •  >ft,+e,-i- 

Note  that  e,-  +  1  is  an  upper  bound  on  /,,  the  nnmber  of  farthest  neighbors  of  p,-. 

Claim  2  The  points Pv,, though  not  necessarily  distinct,  are  in 
sorted  order  relative  to  our  clockwise  ordering  of  the  polygon’s  vertices. 

Proof  Suppose  not,  i.e.  3 i,j  such  that  i  <  j  but  pii  follows  ft,  in  the 
clockwise  ordering  of  vertices.  By  the  quadrangle  inequality,  we  must  have 

<*(Pi»PO  +  d(p„ pi;)  <  d(p,',PbJ)  +  d(p,-,pii). 

But  this  contradicts  our  assumption  that  p*.  is  a  farthest  neighbor  of  p,  and 
ft,  is  a  farthest  neighbor  of  py.  ■ 

Claim  3  For  0  <  i  <  n  -  1,  either  precedes  pfj+1  in  the  clockwise 

ordering  of  vertices  or  is  the  same  point.  In  other  words,  the  ranges  associated 
with  pi  and  Pi+i  may  overlap  only  at  the  point  ft.+I . 

Proof  Suppose  not,  i.e.  3*  such  that  ftj+l  strictly  preceedes  ft,+e,  in  the 
clockwise  ordering  of  vertices.  By  the  quadrangle  inequality,  we  must  have 

rffotft rt-ej  +  ^te+i.ft*,)  <  <*(*•>«,.«) +  <*te+i.  »<+*,)• 

But  this  contradicts  our  assumption  that  p^^  is  a  farthest  neighbor  of  p,- 
and  p^+1  is  a  farthest  neighbor  of  p,+ t.  ■ 


far  the  proof  is  largely  from  the  solution  by  Joel  Weis  to  problem  4  oe  homework 
2  Car  a  cover  on  computational  geometry  (18.40*)  taught  ia  the  spriag  of  1*88  by  Alok  Aggarwal. 
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The  preceding  two  claims  imply  that 


n-l 

E  e«  ^  »• 

t*0 


Thus, 

n— 1 

+  1)  <  2n 

icO 

and 

■ 

n-1 

£  /<  < 2n- 
isO 

Note  that  for  odd  n,  any  regular  n-gon  gives  ns  an  example  of  a  convex  polygon 
for  which  /,  =  2n  since  every  vertex  of  such  an  n-gon  has  exactly  two  farthest 
neighbors. 
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